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ABSTRACT 



An analytical study of linear-amplifying instabilities of a 
laminar boundary layer as found in the experimental data of the 
LaRC/8— Ft. laminar-flow control < LFC) experiment was completed 
and the results are presented. The LFC airfoil used for this 
experiment was a swept, supercritical design which removed suc- 
tion air through spanwise slots. The amplification of small dis- 
turbances by linear processes on a swept surface such as this can 
be due to either Tol lmien-Schl icht ing < TS) and/or crossflow < CF> 
mechanisms. This study consists of the examination of these two 
instabilities by both the commonly used incompressible (SALLY and 
MARIA) analysis and the more involved compressible ( COSAL) analy- 
sis. A wide range of experimental test conditions with varia- 
tions in Mach number, Reynolds number, and suction distributions 
were available for this study. Experimentally determined tran- 
sition locations were found from thin— film techniques and were 
used to correlate the "n-factors" at transition for the range of 
test cases. 



mi" wmmm tin a | HUH i II iHUMi.i II Hi ililIiilUllLJliillllljli jj tiill I n 1 1 , .iIJid.hi.Ls j ilMiHlilil 


BACKGROUND/ MOTIVATION 


The design of laminar-flow airfoils is dependent upon the 
use of transition prediction techniques which accurately model 
the characteristics of the boundary layer. A technique which is 
widely used today is the e 11 method that was developed sepa- 
rately by Smith^ and Van Ingen^. This semi-empirical technique 
examines the tendency of a small disturbance wave within the 
laminar boundary layer to amplify over a given surface to the 
point where breakdown to turbulence occurs. The logarithmic 
amplification ratio, or n-f actor as defined in figure 1, is a 
measure of the spatial growth in amplitude of a disturbance wave 
from the initial point of instability. Stability theory assumes 
that the most amplified disturbance (for a critical combination 
of frequency, wave-angle, and wavelength) is the driving mecha- 
nism behind transition. Based on calibrations of a linear 
incompressible theory with low-disturbance wind tunnel data, 
stability analyses have shown fairly consistently that transition 
corresponds to n-f actors on the order of ten. However, since 
stability theory does not account for the influence of the free- 
stream disturbances, it seems fair to assume that n-factor 
calibrations are only valid between data sets with similar 
disturbance environments^. 

This research was conducted for several reasons: to examine 

the amplification of small disturbances by linear processes on a 
modern swept airfoil; to verify the use of a linear stability 
theory as a design tool through transonic speeds; and to compare 
and contrast the incompressible and compressible methods. 


• Accurate transition prediction techniques are essential for 
laminar-flow airfoil design 

• e n method developed separately by Smith and Van Ingen 

- Semi-empirical technique 

- Relates stability theory (O-S Eq.) to transition 

- Uses amplification ratio as indicator of transition 

• Solutions to 0-S equation give local amplification rates 

n = ln(A/A 0 ) = /* (local amplification rate) dx 
*o 

• Items to note: 

- Solutions depend on proper selection of f, \}/, and X 

- Theory does not account for initial disturbance levels (A 0 ) 

• Objectives of analysis 

- Validate codes used in design (SALLY & MARIA) 

- Examine effect of compressibility (COSAL) 

a 

Figure 1 


473 


BASIC CORRELATION METHOD 


Figure 2 depicts the basic approach employed for this 
research. Measured pressure and suction distributions were used 
to calculate streamwise and crossflow velocity profiles along the 
chord . A stability analysis of these profiles revealed the 
amount of TS and CF disturbance amplification from the initial 
point of instability. Comparisons of the measured transition 
locations with these calculations determined the correlated 
n-factors at transition. 

Solutions to the governing stability eguation, which is the 
celebrated Orr— Sommerf eld (OS) eguation, provide the necessary 
information for analyzing the growth of disturbances within the 
boundary layer. Two of three disturbance variables, freguency 
(f), wavelength (A), or wave-angle < tji > , must be specified with 
the third given by the solution to the OS eguation. Analysis of 
TS disturbances was accomplished using a fixed freguency and 
wave— angle method. For incompressible flow it was assumed that 
the TS disturbances which amplify most are those moving in the 
direction of the free— stream ( ijt = 0). Thus, the incompressible 
TS (SALLY) analysis consisted of analysis of a range of freguen— 
cies with a wave— angle of zero. Since the two-dimensional dis- 
turbance assumption could not be used for compressible calcu- 
lations, the compressible TS ( COSAL) 6 analysis consisted of 
analysis of, first, a range of freguencies with a wave-angle of 
zero, then a range of wave angles for the most critical fre- 
guency. Suffice to say that the COSAL analysis can be guite time 
consuming. The CF analysis was accomplished with the MARIA 7 code 
which gives approximate solutions using a fixed wavelength and 
freguency method. 



iBoundary-lay er. Stability , 

Experimental 1 (CEBECI) (SALLY, MARIA or COSAL)' Transition 

data | analysis | analysis | n-factors 


Figure 2 
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EXPERIMENTAL DATA 


The 8-Foot LFC experiment had all the basic information 
necessary for a stability theory correlation — a pressure and 
suction distribution as well as transition measurements®. 

Because the upper surface was designed as the primary test 
surface, only upper surface data were analyzed. A wide range of 
experimental test conditions with variations in free— stream Mach 
number chord Reynolds number <R ), and surface suction < C a > 

were available for this study. As shown in figure 3, the avail- 
able data encompassed both the subsonic and transonic regimes. 

For the subsonic cases, the pressure distributions all had a 
strong leading-edge peak as is shown in figure 3 by a represen- 
tative case at M w = 7. The suction distributions available for 
these cases were all of the type designed for the airfoil with 
suction over the full chord. Only the overall suction levels 
were changed between cases without any attempt to optimize the 
distributions. At moderate to high overall suction levels, full- 
chord laminar flow was achieved for the entire subsonic range of 
H k and R c> 

The transonic cases available were all at the design Mach 
number < M k = .82). As is shown in figure 3, the pressure 
distributions for these cases closely followed the design distri- 
bution except for the slight waviness and slightly higher velo- 
cities towards the nose. The waviness in the pressure distri- 
butions was thought to be due in part to the model deflecting 
under loads. A sonic bubble encompassed most of the upper 
surface and culminated in a weak shock at the higher Reynolds 

numbers < R„ > 12 million). Full-chord laminar flow was achieved 
c 

until the shock appeared and then transition occurred ahead of the 
shock. Because of this, analysis was done on the "hybr id"-type 
suction distributions in which suction was systematically shut 
off at the trailing edge, moving towards the leading edge. 


Upper surface 


SUBSONIC TRANSONIC 

. 4 , . 6. .7 .82 
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SAMPLE SUBSONIC CALCULATIONS 


Figure 4 shows the results of an incompressible TS analysis 
of a sample subsonic case ( = .6 and R c = 10 million). For the 
pressure and suction distributions shown in the figure* n— factors 
were calculated for a range of frequencies with a wave-angle of 
zero. The combination of the strong adverse pressure gradient 
near the leading edge and the full-chord suction distribution 
confined the TS amplification to within the first 15% chord 
followed by a stable region to the trailing edge. For this case* 
with its relatively low overall suction level* transition was 
measured to be around 10—15% chord. As is shown in figure 4 , 
the TS growths reached a logarithmic amplification ratio of 
n = 10 in this region, for the cr it ical frequency of f = 10 kHz. 


Mqq- .6, R c = 10 million 
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Figure 4 



SAMPLE SUBSONIC CALCULATIONS 


Figure 5 shows the results of an incompressible CF analysis 
of the sample subsonic case. The analysis indicates that CF 
instability is not strong enough to cause transition in this 
case. In fact, over the entire subsonic range of M w and R 
cases available, CF logarithmic amplification ratios did not c 
exceed n = 2.5 in the nose region where the cross-flow instability 
is considered to be most critical. 


Moo”. 6, Rc = 10 million 
INCOMPRESSIBLE CROSS FLOW (MARIA) 



Chord station, x/c 


Figure 5 
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RESULTS OF SUBSONIC ANALYSIS 


The results of the subsonic incompressible analysis are 
summarized in figure 6 which shows the calculated transition 
n-factors versus the subsonic range of Mach numbers for the 
Reynolds numbers of 10 and 20 million. Since CF amplification 
was shown to be ins ignif icant f transition was due entirely to TS 
amplification. The figure shows the TS logarithmic amplification 
ratios that correlated with transition. The TS transition n— factors 
varied from n = 8. 5 to 10.5 over the range of M k and R c ana- 
lyzed. This variation in n-factor is attributed to the fact that 
the wind tunnel disturbance environment is changing with M^ and 


INCOMPRESSIBLE TS 



Figure 6 
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CONCLUSIONS OF SUBSONIC ANALYSIS 


The main conclusions that were drawn from the incompressible 
subsonic analysis 9 are listed in figure 7. The transition pro- 
cess was found to be dominated by TS amplification. For these 
subsonic cases, transition occurred when TS logarithmic amplifi- 
cation ratios in the range of n = 8. 5 to 10.5 was reached near the 
leading edge. 


• CF disturbance amplification at leading edge 

- Small 

- Unaffected by suction 

• Decreasing overall suction levels 

- Decreased amount of laminar flow 

- Allowed higher TS growths at leading edge 

• Therefore, transition process 

- Dominated by TS disturbances 

- Occurred when njs in the range 8.5 to 10.5 
was reached 


Figure 7 
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SAMPLE TRANSONIC CALCULATIONS 


Figure 8 shows the results of an incompressible TS analysis 
of a sample transonic case < M w = .82 and R c = 20 million). For 
the pressure and suction distribution shown in the figure* 
n-factors were calculated for a range of frequencies with a wave- 
angle of zero. The combination of the wavy pressure distribution 
and suction being turned off after 8% chord allowed strong TS 
amplification to begin. For this "hybrid" case, transition was 
measured to occur between 20 and 28% chord. As is shown in the 
figure, TS disturbances grew to a logarithmic amplification ratio 
of n = 10 to 13 in this region, for a critical frequency of 
f = 9 kHz. 
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Figure 8 
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SAMPLE TRANSONIC CALCULATIONS 


Figure 9 shows the results of an incompressible CF analysis 
of the sample transonic case. As is shown in the figure* CF 
disturbances grew quickly to a peak logarithmic amplification 
ratio of n = 4. 5 at about 5% cord and then rapidly damped out. 
There was no calculated CF amplification in the region of 
measured transition. For this reason* it was concluded that CF 
disturbances did not contribute to transition. The case shown 
here illustrates the "worst" case calculations in terms of CF 
disturbances over the entire range of M c _ and R . 


M c0 = .82 R c =20 million 

INCOMPRESSIBLE CROSS FLOW (MARIA) 
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VARIATION OF TRANSITION FOR SIMULATED HLFC 


During the experiment, an attempt was made to simulate 
possible "hybrid" laminar-flow control (HLFC) conditions by 
turning off the suction level at various chord locations. The 
HLFC cases consisted of measuring the extent of laminar flow as 
the upper surface suction was systematically shut off, at the 
same time maintaining lower surface suction on only the first 25% 
chord. This simulated HLFC approach is considered somewhat "non- 
ideal istic" since no effort was made to seal or smooth the un- 
sucked regions of the airfoil. Furthermore, it is clear that the 
current pressure distribution (measured or theory) is not like 
that which would be designed for HLFC. However, these cases 
provide data for which transition occurs far forward of shock 
location thereby being ideal for a stability theory correlation 
at transonic speeds. 

Figure 10 shows the measured variation of transition 
location on the LFC airfoil upper surface with extent of suction 
for M k = .82 and two chord Reynolds numbers. For R c = 10 
million, transition occurred roughly 20% chord past the point 
where suction was shut off, until full-chord laminar flow was 
achieved. For R = 20 million, transition occurred roughly 10% 
chord past the point where suction was shut off, until the shock 
present on the rearward end of the model tripped transition. The 
data clearly indicate the extent of laminar flow that could be 
maintained beyond the point where suction stopped for transonic 
airfoils with weak— shock and shockless conditions. 


SWEPT LFC AIRFOIL SLOTTED UPPER SURFACE 


L.S. suction, O^x/ci.25 



Figure 10 



TRANSONIC RESULTS 


Figure 11 indicates the results of the incompressible 
stability analysis of the HLFC data that has shockless flow on 
the upper surface < R^ = 10 million). The figure shows the range 
of TS amplification at transition versus the extent of suction 
applied to the airfoil. For each case in which the suction was 
turned off at a different chord location, symbols were plotted 
for the calculated incompressible TS logarithmic amplification 
ratio that corresponded with the measured transition zone. The 
open symbols indicate the n — factor calculation that corresponded 
to the chord location of the last fully laminar thin- f ilm gage, 
while the closed symbols indicate the calculation that corre - 
sponded to the first fully turbulent gage. The transition zone 
shown on the figure helps illustrate the correlated n-factors 
for the extent of suction used. For this model at test condi- 
tions M w = .82 and R c = 10 million, the TS transition 
n-factors for cases with suction past 30% chord are slightly 
lower than those found with little or no suction. The cases with 
the extent of suction less than 30% chord have incompressible TS 
amplification ratios at transition of n = 10. The cases with 
suction past 30% chord have correlated n-factors closer to n = 7. 
P®rhaps this indicates that the transition process on swept 
LFC/HLFC airfoils becomes increasingly non-linear as the extent 
of desired laminar izat ion increases. 


Calculated n-factors lor variable chordwise suction 
M„ .82 R c - 10 X 10 6 



Figure 11 
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TRANSONIC RESULTS 


Figure 12 indicates the results of the incompressible 
stability analysis of the HLFC data that has a weak shock formed 
on the rear of the upper surface <Rc = 20 million). As is shown 
in figure 11 the range of TS amplification at transition is 
plotted against the extent of suction applied to the airfoil. 

the cases ranging from no suction to suction over the first 
30% chord, the same trend of correlated n-factors of n = 10 was 
found. Suction beyond 30% chord, however, shows that transition 
is increasingly being dominated by non-linear processes as the 
extent of suction approaches the shock location. These results 
demonstrate that for a wind tunnel environment, transition can be 
very sensitive to the formation of a shock. 


\ 


Calculated n-factors for variable chordwise suction 
Moo = 82 R q = 20 X 10 6 
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Figure 12 
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CONCLUSIONS OF TRANSONIC ANALYSIS 





/ 




: ' The conclusions that were drawn from the incompressible 

analysis of the transonic data are listed in figure 13. For 

cases in which full-chord suction distributions were used, 

transition was sensitive to the rearward attachment point of the 

sonic bubble. If shockless flow was attained < R_ < 12 million) 

c 

then transition could be maintained very close to the trailing 
edge. However, if a weak shock appeared < R c > 12 million) then 
transition occurred ahead of the shock for a wide range of over- 
all suction levels which could not be explained with a linear 
stability theory. Turning off suction ahead of shock, which 
simulated an HLFC concept, allowed transition to occur as a conse- 
quence of disturbance wave amplification and good theory correla- 
tions were found for the available range of R c> The CF disturbances 
were generally small and unaffected by suction, while TS disturb- 
ances grew rapidly to transition once suction was turned off. 

The incompressible n-factor correlations at transition for this 
transonic data were n = 10 when suction was applied up to 30% 
chord and n = 7 when suction was applied further than 30% chord. 


• Full-chord suction cases 

- Transition occurs ahead of shock location for R C M2 million 
which could not be explained by linear stability theory 

- Transition occurs at trailing edge for shockless flow (R c < 12 million) 
•"Hybrid" suction cases 

- Turning off suction ahead of shock allowed natural transition and 
good theory correlations for range of R c 

- CF disturbances small and unaffected by suction 

- TS disturbances grew rapidly to transition after suction turned off 

• Incompressible correlations at transition 

- 10 near leading edge 

- 7 in the mid-chord region 


Figure 13 



485 


SAMPLE TRANSONIC CALCULATIONS 



Figure 14 shows the results of a compressible TS analysis of 
the sample transonic case (M tt = .82 and R c = 20 million). For 
the pressure and suction distribution shown, compressible 
n-f actors were calculated first for a range of frequencies with a 
wave angle of zero, then for a range of wave angles for the most 
critical frequency. Since transition was measured between 20 and 
28% chord, the compressible TS n-f actor at transitions corre- 
sponds to n = 5 - 7.5 for the most critical disturbance 
( f = 7 kHz and \Jj = -50°). 
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COMPARISONS OF SALLY AND COSAL RESULTS AT TRANSITION 


Comparisons of incompressible and compressible stability 
analysis methods were conducted for a range of M and the 
results are tabulated in figure 15. As shown in the figure, 
compressibility effects reduced the calculated n-f actors at 
transition. The difference accounts for nearly 25% reduction of 
t he calculated amplification for t he — 0.4 to over 50% 

reduction at M w = .82. Both have fairly consistent correlations 
with the incompressible n-f actors of 9 to 10 and the compressible 
n-f actors of 5 to 6. However, as noted earlier, the COSAL anal- 
ysis involves solutions based on a second unknown <^> which 
increases the time and expense of each compressible case. 
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• Both have fairly consistent correlation 

• COSAL introduces another unknown - M 

• COSAL is more expensive and time consuming 


Figure 15 
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SUMMARY 


A wide range of low-disturbance wind tunnel data with 
variations in M^, R c and suction were analyzed based on a 
linear boundary-layer stability theory. The data, which were from 
a large chord, swept, supercritical LFC experiment, encompassed 
both the subsonic and transonic flow regimes. The analysis 
showed that as far as 1 aminar izat ion through the control of dis- 
turbance amplification, the model performed as designed. The CF 
amplification was successfully controlled through design by using 
a proper combination of Reynolds number, pressure gradients and 
sweep angle. The TS amplification could be successfully controlled 
by suction as long as strong non-linearities <e.g. shocks) were 
not present on the airfoil. As is shown in figure 16, the TS 
logarithmic amplification ratios that were calculated at transi- 
tion were n = 9 to 11 for incompressible theory and n = 4 to 6 
for compressible theory. Based on these results, it seems fair 
to conclude that linear stability theory can be used as a guide 
in the design of laminar-flow airfoils through transonic speeds. 
However, a word of caution is needed, increasing conservatism in 
the choice of a design n— factor is required with increasing 
lengths of LFC, especially at transonic speeds. 


• For the range of subsonic and transonic cases 
TS disturbances dominated the transition process 

• Transition corresponded to 

- Incompressible values of 9^nyg^11 

- Compressible values of 5^nyg^6 

• Linear stability theory can be used as a guide in 

the design of laminar flow airfoils at transonic speeds 

• Increasing conservatism in choice of n-factor is required 
with increasing length of LFC 


Figure 16 
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